Valley contrasting physics in graphene: magnetic moment and topological transport 
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We investigate physical properties that can be used to distinguish the valley degree of freedom in 
systems where inversion symmetry is broken, using graphene systems as examples. We show that 
the pseudospin associated with the valley index of carriers has an intrinsic magnetic moment, in 
close analogy with the Bohr magneton for the electron spin. There is also a valley dependent Berry 
phase effect that can result in a valley contrasting Hall transport, with carriers in different valleys 
turning into opposite directions transverse to an in-plane electric field. These effects can be used to 
generate and detect valley polarization by magnetic and electric means, forming the basis for the 
so-called valley-tronics applications. 

PACS numbers: 73.63.-b,75.75.+a,85.35.-p 



Graphene, the monolayer carbon honeycomb lattice, 
has extraordinary electronic properties [TJ [5J. Its band 
structure has two degenerate and inequivalent valleys at 
the corners of the Brillouin zone. Because of their large 
separation in momentum space, inter- valley scattering is 
strongly suppressed O 0] , implying the potential use of 
valley index in a way similar to the role of spin in spin- 
tronics applications. Interesting valley dependent phe- 
nomena are being actively explored [5]. 

In this Letter, we propose a general scheme to gener- 
ate and detect valley polarization in graphene systems 
with broken inversion symmetry. We reveal that there 
is an intrinsic magnetic moment associated to the valley 
index, in close analogy with the Bohr magneton to the 
electron spin. This property makes the valley polariza- 
tion a directly measurable physical quantity. The broken 
inversion symmetry also allows a valley Hall effect, where 
carriers in different valleys flow to opposite transverse 
edges when an in-plane electric field is applied. It opens 
a new possibility to the much desired electric generation 
and detection of valley polarization. The valley Hall ef- 
fect is analogous to the spin Hall effect [B], and falls into 
the same category as the Berry-phase supported topolog- 
ical transport phenomena. 

Graphene systems with broken inversion symmetry are 
of direct experimental relevance. Zhou et al. [7] have re- 
cently reported the observation of a band gap opening in 
epitaxial graphene, attributed to the inversion symmetry 
breaking by the substrate potential. In addition, in bi- 
ased graphene bilayer, inversion symmetry can be explic- 
itly broken by the applied interlayer voltage More- 
over, as we show below, the emergent valley contrasting 
physics is a generic consequence of bulk symmetry prop- 
erties, which provides a new and much standard pathway 
to potential applications of ' valley tronics' in a broad class 
of semiconductors [TO], as compared to the novel valley 
device relying on the peculiar property of the edge state 
in graphene nanoribbon [5 . Graphene with broken in- 
version symmetry serves as a paradigm to demonstrate 
the general features and necessary conditions of such ap- 



plications. 

Before starting specific calculations, it will be instruc- 
tive to make some general symmetry analysis. A valley 
contrasting magnetic moment has the relation m„ = x T z , 
where t z — ±1 labels the two valleys and x is a coeffi- 
cient characterizing the material. Under time reversal, 
m,, changes sign, and so does t z (the two valleys switch 
when the crystal momentum changes sign). Therefore, 
X can be non-zero even if the system is non-magnetic. 
Under spatial inversion, only r 2 changes sign. Therefore 
m„ can be nonzero only in systems with broken inversion 
symmetry. 

Inversion symmetry breaking simultaneously allows a 
valley Hall effect, with j v — u v H z x E, where a v H is the 
transport coefficient (valley Hall conductivity), and the 
valley current j v is defined as the average of the valley 
index times the velocity operator. Under time reversal, 
both the valley current and electric field are invariant. 
Under spatial inversion, the valley current is still invari- 
ant but the electric field changes sign. Therefore, the val- 
ley Hall conductivity can be non-zero when the inversion 
symmetry is broken, even if the time reversal symmetry 
remains. 

Armed with the insight from the above symmetry 
analysis, we now consider a concrete example, a sin- 
gle graphene layer with a staggered sublattice potential 
breaking the inversion symmetry. Staggered sublattice 
potential is generally expected in epitaxial graphene as 
pointed out in the review by Geim and Novoselov [2] 
and explicitly shown by ab initio studies In the 

tight binding approximation, it can be modeled with a 
nearest-neighbor hopping energy t and a site energy dif- 
ference A between sublattices [T2], [T3J [T^]. For relatively 
low doping, we can resort to the low-energy description 
near the Dirac points. The Hamiltonian is given by 

H = — at (q.x T z°~x + q y o- y ) + —a z , (1) 

where er is the Pauli matrix accounting for the sublat- 
tice index, and q is measured from the valley center 
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FIG. 1: (color online). Energy bands (top panel) and orbital 
magnetic moment of the conduction bands (bottom panel) of 
a graphene sheet with broken inversion symmetry. The Berry 
curvature fi(fe) has a distribution similar to that of m(fc). 
The first Brillouin zone is outlined by the dashed lines, and 
two inequivalent valleys are labeled as Ki and K2. The top 
panel shows the conduction (red) and valence (blue) bands in 
the energy range from —1 to 1 eV. The parameters used are 
t = 2.82 eV and A = 0.28 eV. 



Ki,2 = (t47t /'Sa)x with a being the lattice constant. 
In the following we shall focus on the n-doped graphene. 
Generalization to the p-doped graphene is straightfor- 
ward due to the particle-hole symmetry presented in this 
system. 

Because spin-orbit coupling is extremely weak in 
graphene [TS], the valley magnetic moment can only be 
of orbital nature. To study this quantity, we invoke 
the semiclassical formulation of the wavepacket dynam- 
ics of Bloch electrons [16]. It has been shown that in 
addition to the spin magnetic moment, Bloch electrons 
carry an orbital magnetic moment given by tn(fe) = 
-i(e/2h)(Vku\ x [H(k) - e(fc)]|V fe zt), where |u(fc)) is 
the periodic part of the Bloch function, H(k) is the Bloch 
Hamiltonian, and e(k) is the band energy |16j . It orig- 
inates from the self-rotation of the wavepacket. For a 
two-dimensional system, the orbital magnetic moment is 
always in the normal direction of the plane and may be 
written as m(k)z. Its momentum dependence can easily 
be calculated from the tight-binding Bloch states, and is 
shown in Fig. [I] As we can see, rrt(fc) is concentrated 
in the valleys and has opposite signs in the two inequiv- 
alent valleys. Analytic expression can also be obtained 
from the model Hamiltonian (1) in the neighborhood of 
such valleys: 
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where m* = (2Ah 2 )/(3a 2 t 2 ) is the effective mass at the 
band bottom. This is in close analogy with the Bohr 
magneton for the electron spin, where the effective mass 
becomes the free electron mass. In fact, the analogy goes 
further, because one can also obtain the spin Bohr mag- 
neton by constructing a wavepacket at the bottom of the 
positive energy bands of the Dirac theory and calculat- 
ing the self-rotating orbital moment. Therefore, it makes 
sense to call the orbital moment calculated above as the 
intrinsic magnetic moment associated with the valley de- 
gree of freedom, provided one is only concerned with low 
energy electrons near the bottom of the valleys [T51 [T5] . 

The valley magnetic moment has important implica- 
tions in valleytronics as it can be inferred from all kinds 
of experiments analogous to those on the spin magnetic 
moment. For example, while spin polarization of elec- 
trons can be created by a magnetic field (Pauli para- 
magnetism), we expect a similar valley polarization in 
graphene due to coupling between a perpendicular mag- 
netic field and the valley magnetic moment. Moreover, 
for typical values of A ~ 0.28 eV and t - 2.82 eV with 
a lattice constant a = 2.46 A we find (i* B to be about 30 
times of the Bohr magneton. Therefore the response to a 
perpendicular magnetic field is in fact dominated by the 
valley magnetic moment at low doping in graphene. In- 
terestingly, unlike the spin moment which will respond to 
magnetic fields in all directions, n* B only couples to mag- 
netic fields in the z-direction. This strong anisotropic 
magnetic response may be used to distinguish the spin 
and valley magnetic moment. 

Complimentarily, a population difference in the two 
valleys may be detected as a signal of orbital magneti- 
zation. The orbital magnetization consists of the orbital 
moments of carriers plus a correction from the Berry cur- 
vature l20l 
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[m(k) + (e/h)(fx-s(k))Q(k)], (4) 



It is instructive to consider the low energy limit (q 



where fi is the local chemical potential, and the inte- 
gration is over states below the chemical potential. The 
Berry curvature fi(fe) = Q,(k)z is defined by fi(fe) = 
Vfc x (it(fe)|£V k\u(k)) and its distribution has a sim- 
ilar structure to that of m(fe). We note that Eq. Q 
is for temperatures much lower than the energy scale 
of band structure (roughly given by A), which holds 
up to room temperature as the experimentally observed 
bandgap A ~ 0.28 eV [J\. For two-band model with 
particle-hole symmetry, we have a simple relation be- 
tween the orbital magnetic moment and the Berry cur- 
vature in the conduction band: m(fc) = (e/H)e(k)Cl(k). 
Using this relation, Eq. Q may be further simplified as 
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FIG. 2: (color online). Electric generation (a) and detection 
(b) of the valley polarization, (a) An in-plane electric field will 
generate a transverse valley current, which leads to a net val- 
ley polarization on the sample edges, (b) A valley-polarization 
created by the valley filter [5] results in a transverse voltage 
across the sample. 



M = 2(e/K) J J4j^0(fe). When the two valleys are 
in equilibrium (the chemical potential fi is common to 
both), this integral vanishes because the Berry curvature 
has opposite values in the two valleys. In the presence of 
a population difference, the chemical potential has differ- 
ent values in the two valleys ^ fi 2 - Therefore, the net 
orbital magnetization is given by 



SM = 2-[^iCi(/ii) + ^2C 2 (/U 2 ) 



2-d(fi)5[i , (5) 



where 2nCi(fi) = / Ml d 2 kQ(k) is the Berry phase around 
the Fermi circle in valley Kj, Sfi = [i\ — /i 2 and 2/2 = 
A'i + fJ-2- The approximate equality holds for \n > A, 
where the Berry phases approach ±n. Thus in a crude 
estimation, SM reduces to (e/h)Sfi. 

Next we discuss the Berry-phase supported topological 
transport in our system. It has been well established that 
in the presence of an in-plane electric field, an electron 
will acquire an anomalous velocity proportional to the 
Berry curvature in the transverse direction |16j . giving 
rise to an intrinsic contribution to the Hall conductiv- 
ity < = 2(e 2 A) I $*f(k)n(k), where f(k) 

is the Fermi-Dirac distribution function, and the factor 
of 2 comes from spin degeneracy. There is also a side- 
jump contribution [23 proportional to the Berry curva- 
ture when carriers scatter off an impurity potential. The 
aforementioned symmetry argument manifests itself in 
the symmetry property of the Berry curvature r2(fe): it 
is an odd function in the presence of time reversal sym- 
metry and even in the presence of inversion symmetry. 
From Eq. (JT|) we have for the conduction band 
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Ignoring skew-scattering and other effects due to inter- 
valley scattering, we find a valley-dependent Hall con- 
ductivity as 
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where qp is the Fermi wave vector which is related to 
the bulk chemical potential by /i = ^-\/A 2 + 3q F a 2 t 2 . 
The third term is the side-jump contribution, which is 
also independent of the scattering rate [53]. Interestingly, 
when the Fermi energy Ep — /i is bigger than the gap 
A, such that the Berry curvature peak is well covered 
by occupied states, the Hall conductance approaches a 
quantized value of r z e 2 jh. 

The valley dependence in the Hall current will lead 
to an accumulation of electrons on opposite sides of the 
sample with opposite valley index (see Fig. [2^l). If an 
electric field E y is applied along a strip of the sample, 
the valley population difference at one edge is given by 

Sn = f x T v = <j v H E y T v , <r v H = T z o-H(T z )/e , (8) 

where t v is the inter-valley life time. The valley polar- 
ization is distributed along the edge within the diffusion 
length lp = vp where up is the Fermi velocity 

and To is the intra- valley scattering time. From Ref. [4j 
we take r = 0.1 ps and t v — 50 ps. Assuming an electric 
field E — 1 mV//im, we find a valley population differ- 
ence of 10 — 100 per /j,m along the edge and distributed 
over a width of lp ~ 1/im. This valley polarization may 
be detected as a magnetic signal as we discussed before. 

Clearly, if there is a net valley polarization (/zi ^ 1x2), 
a Hall current will appear upon the application of an 
electric field E,,, 
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This Hall current will then lead to a measurable trans- 
verse voltage across the sample. If the width of the bulk 
region is smaller or comparable to the mean free path, the 
transverse voltage along the edge gives a local mapping 
of the valley polarization in the bulk. We show in Fig. [2)3 
an experimental setup in conjunction with the valley fil- 
ter device [5] to demonstrate this effect (we note that 
inversion symmetry breaking does not change the edge 
state property needed for the valley filter to function). 

The valley magnetic moment and valley Hall effect pre- 
dicted above are generic features in systems with broken 
inversion symmetry, as shown by another example, the 
biased bilayer graphenc. This system may be modeled 
by an intra-layer nearest neighbor hopping t, an inter- 
layer nearest neighbor hopping tj_ , and an energy bias A 
between the layers, which breaks the inversion symme- 
try [5] . Angle-resolved photoemission spectroscopy stud- 
ies [8 of bilayer graphene films synthesized on SiC sub- 
strates confirm the band structure from this model. 
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FIG. 3: (color online). Valley contrasting properties of con- 
duction bands in biased graphene bilayer. (a) Energy dis- 
persion (green or gray curves) and Berry curvature (black 
curves), (b) Orbital magnetic moment. Solid curves for 
lower conduction band and dashed curve for upper conduction 
band. The quantities are shown for the Ki valley. Distribu- 
tions of m(fc) and Q(k) have opposite signs in the K2 valley. 
The corresponding valley magnetization (c) and the valley 
Hall conductivity (d) are also shown as a function of chemical 
potential. The parameters used are t — 2.82 eV, A = 0.2 eV, 
and t± = 0.4 eV. 



Biased bilayer graphene has two positive energy bands 
(conduction) and two negative energy bands (valence) if 
spin degeneracy is discounted. In Fig. [3] we show nu- 
merically calculated energy bands, Berry curvatures and 
orbital magnetic moments of the two conduction bands. 
The parameter values are chosen in accordance with ex- 
perimental result [5]. fi(fe) and m(fc) are again peaked 
at the valley bottom. The valley magnetization and the 
valley Hall effect are of the same order of magnitude as 
in the epitaxial single-layer graphene. We note that the 
valley-dependent Hall conductance approaches a quan- 
tized value of 2r 2 e 2 /h, twice of that for the single layer. 
This is consistent with the fact that in bilayer graphene, 
the Berry phase acquired by an electron during one circle 
around the valley becomes ±2ir instead of ±7r when the 
gap closes [25] . 
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